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^ Abstract 

^ Riemannian cubics are critical points for the norm of acceleration of curves in 

^ Riemannian manifolds M. In the present paper the L°° norm replaces the norm, 

^ and a less direct argument is used to derive necessary conditions analogous to those 

^ in [T7] . The necessary conditions are examined when M is a sphere or a bi- invariant 

^ Lie group. 
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1. Introduction 

j> Consider the task of moving a physical object between given configurations at 

^ given times, in a way that is optimal in some sense yet to be defined. 

^ Example 1. Let t^ < ti be given real numbers, and let Xq,Xi G be given, where 

^ denotes Euclidean m-space. The task is to move a point of unit mass from Xq to 

O Xi along a curve x : [toj^i] — ^ E^ so as to minimise the total kinetic energy 



ti 
to 



K{x) := ^ / \\x^'\t)\\l dt 



X 

where x^^^ denotes the first derivative of x (the velocity), and \\ \\e is the Euclidean 
norm. There is a unique minimising path, namely the uniform-speed interpolating 
line segment. 

Therefore, given real numbers to < ti < . . . < tn and points xo,xi, . . . ,Xn G E^, 
the unique minimal energy interpolating curve is a track-sum of n line segments, 
traversed at uniform speed within each tj\ where j = 1,2, . . . ,n. 
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In practice there may he problems if an object is moved in this way, due to instan- 
taneous changes of momentum at junctions tj for j ^ 0,n. So another measure of 
optimality might be used instead of K . A simple alternative is the mean squared 
norm of the acceleration 

Ux) := / ||x(2)(t)f dt 

ti - to Jta 

where x varies over piecewise-C^ interpolating curves. As seen in 18], the natural 
cubic spline is the unique minimiser of J2. 

For objects more complicated than a point mass, similar considerations may hold. 

Example 2. Suppose that we want to move a rigid body whose centre of mass is 
fixed through given configurations at times tQ,ti, . . . ,tn. The trajectory becomes a 
curve X : [tc^n] M, where M is the 3-dimensional rotation group S0{3). The 
given configurations are points Xj G 50(3), and K is now given by 

K{x) = I f'Wx^^Wl^^dt 

J to 

where \\ \\z is the Riemannian norm at z & M . As in the Euclidean case, a minimiser 
of K is a track-sum of geodesic arcs traversed at uniform speed within each [tj^i,tj]. 

Such a path might be impractical, for the same reasons as in Example [I| So, as 
before, we consider minimisers of 

J2{x) := / ||Via;(i)||^(,) dt 

where VtX^^^ is the covariant acceleration of the interpolant x : [to;^n] M, and is 
adapted to the Riemannian structure of M . Minimisers of J2, called Rieman cubic 
splines, are track-sums 0/ Riemannian cubics defined over the tj-i,tj, with trivial 
covariant accelerations at to one? t„. 

Riemannian cubics are curves that with more mathematical structure than either 
geodesies or cubic polynomials, and are given by the nonlinear ordinary differential 
equation |77| / 

V,3x« + i?(Vtx(i\a;(i))a;(^) = (1) 
where R is the Riemannian curvature tensor field. 
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Riemannian cubics are reviewed in [23] • More details, including recent results, can 
be found in [18], [19], [20], [21], [22], [21], [26], [H], [6], [7], [12]. The present paper 
initiates an alternative line of enquiry, motivated by different engineering considera- 
tions. 

Whereas J2 measures the mean squared norm of the force applied to the moving 
object, it may sometimes be more important to minimise the maximum norm of 
applied force, namely 

Joo(x) := max{||Vtx(^)||a;(i) : t G [to,^n]} 



where x : [to, tn] ~^ ^ varies over piecewise-C^ interpolants. As seen in ^10, in order 
to find minimising curves it suffices to take n = 1 and impose additional constraints 
on velocities at endpoints. 

More precisely, given to < ti together with G M and vi G TM^^, a piecewise-C^ 
curve X : [to, ^1 

is FxQ^xi,v(j,vi'f(^(^sihle when x{ti) — Xi and x^^^{ti) — Vi for 
i = 0,1. A minimiser of Joo on Fa^g^^^j^^jg^^^ -feasible curves is said to be FxQ^xi,vo,Vi_- 
optimal. Our main result. Theorem [T| gives a necessary condition for a curve to 
be FxQ^xi,vo,vi-optimaA. Examples are given in ^and ^for the cases where M is a 
sphere and a bi-invariant Lie group respectively. 



j9] introduces the related notion of Fxg^xi,vo-optimality, for which an analogous nec- 
essary condition is given by Theorem [2| Then Corollary |2] of ^ 10 says that, if x 



minimises Joo over piecewise-C^ interpolants of many points, then a least one track- 
summand of x satisfies the conclusions of Theorem [T] or Theorem [2l 



2. Optimal Curves 

Let M be a C°° manifold of dimension m > 1, with a C°° Riemannian metric 
( , ). Let real numbers to < ^ibe given. The derivative x^^^ of a piecewise-C^ curve 
X : [to,ti] — )■ M is piecewise-C^, and the Levi-Civita covariant acceleration VtX^^^ is 
a piecewise-continuous field defined along x. 

Definition 1. A piecewise-C"^ curve x : [to,ti] M is said to he feasible when 
x(ti) = Xi and x^^^{ti) = vi for z = 0, 1. 

Let F 

— FxQ^xi,VQ,vi — FxQ^xi,vo,vi,to,ti be the set of all feasible curves, and define 
J : F M by J{x) := max{|| Vta;(^^||a;(t) : t G [to,ti]} where, for w G M, || ||^ is the 
Riemannian norm on TM^ associated with the Riemannian inner product { , )w ■ 
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A minimizer of J is said to be optimal. Our main result is a necessary condition for 
optimality, in terms of 

L{X) := Vi'X + i?(X,x(i))x« 
where R denotes Riemannian curvature, and X a piecewise-C^ field defined along x. 

Theorem 1. Let x E F be optimal. Then \\VtX^^^\\x{t) is constant. For some 
piecewise-C^ function $ : [^07^1] — ^ [0, c>o) with $~-^(0) discrete, we have 

L{<pWtX^^^) = (2) 

where (j){t) := a/ ^{t)- Moreover x and are C°° except possibly where ^(t) = 0. 

Example 3. Let x be an optimal curve in an open convex subset M of Euclidean 
m-space E"^. Then plj reads 



= ^ = A + Bt 

dt'' 

where A, B E E"^ . We have z(j){t) = \\A + Bt\\ where z is the constant length of x^'^\ 

• If z = then A = B = and x is an affine line segment. 

• If A -\- Bt2 = for some t2 E [to, ti] where B ^ then 

^^^^ ^ \t-t2\\\B\\ ^ ^(2) ^ 2: sign(t-t2)^ fort^t2. 
a \\B\\ 

So X is a quadratic polynomial spline which is C°° except possibly at ^2- 

• If A and B are linearly independent then is never zero. By shifting [to;^i] we 
can suppose {A,B) = 0. Setting a := and /3 := ||-B||, 

x^'\t) = z- ^^^^ 



^0? + /32t2 



xit) = z ^ — ^ A + 

a' log(/3t + y/g^ + ^H^) + ^t^a^ + ^ , , ^ 
z — B + Ct + D 

where C,D E E"^ are constants. 



The conclusions of Example [3] also follow from the argument preceding the proof 
of Theorem 4.1 in §4 of [I5]. However [15], uses linear methods from classical ap- 
proximation theory, building on [9], [lO], [H] among others. These methods are 
inapplicable in a more general Riemannian setting. For instance, in ^ where M is 
a sphere, or in ^ where M is a bi-invariant semisimple Lie group, there seems to be 
no alternative to Theorem [H 

The proof of Theorem [T] is less direct than the derivations of the analogous Euler- 
Lagrange equations for geodesies elastic curves [27], and Riemannian cubics 
[T7] . As with these other variational problems, we proceed by considering extremals, 
where an optimal curve is necessarily extremal, but we need two distinct notions of 
extremal. 

In ^ so-called L°°- extremals are used to prove that || VtX*^"'^''(t)||2;(t) is constant, and 
that optimal curves are locally optimal. After some preliminaries in ^ local optimal- 
ity is used in ^ to reduce Theorem [T] to the special case where M is an open subset 
of M™ with a nonstandard Riemannian metric. For this special case, constancy of 
II Vja;''"'^-'(t)||2;(t) is used to identify x with an extremal in the sense of optimal control 
for a control system with state variables in M^™"*"^ and controls in the Euclidean unit 
sphere 5*™"^. Then in ^ a local calculation based on the Pontryagin Maximum 
Principle [2S] is used to prove Theorem [l] for the special case, and this completes the 
proof of Theorem [T] in general. 

A more direct proof of Theorem [Tj is possible for optimal curves whose image is 
contained in a single coordinate chart, based on the Pontryagin Principle for a slightly 
different optimal control problem where ||VtX*^^'*||2,.(t) is not assumed to be constant. 
Then ^ could be omitted, at the expense of generality. Another alternative might be 
to use the Pontryagin Principle for multiprocesses [5], but [|3] seems less complicated. 

It might be thought that an application of the geometric formulation of the Pontrya- 
gin Principle could replace the local calculations of ^ and perhaps do away with 
[|3} For this to work, the space of controls would need to vary from point to point, 
which is not envisaged in [3], [2], [1]. Possibly this might be done in the framework 
of [1] , but our local calculation has the advantage of being more elementary. 

Having proved Theorem [T| we go on to investigate some consequences. In ^ M is 
taken to be the Euclidean unit m-sphere, and the differential equation ^ is rewritten 
as the system Q, (16), without covariant derivatives. Figure [l] in Example |4] shows 



a numerical solution for x. 

In [|8} M is taken to be a Lie group G with a bi-invariant Riemannian metric. The 
differential equation ^ reduces and integrates to a differential equation (18) in the 
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Lie algebra Q, in terms of a Lie reductions V,Xl. Another two conserved scalar 
quantities in terms of V, Xl are also noted, and it is shown that generically x can be 
recovered by quadrature from V,Xl. A class of so-called null optimal curves turn 
out to be the same as null Riemannian cubics [18], [20], [22] and Chapter 4 of [Mj . 
The curious geometry of reductions of non-null optimal curves in 5*0(3) is illustrated 
by Figures [2| [3| |4] of Example [5} 

In ^ M is once more a general Riemannian manifold, and the notion of feasibility 
is relaxed, so that x^^^ is specified only at one endpoint. Simple modifications to the 
proof of Theorem [T] result in an additional necessary condition, given in Theorem [2] 

In [ p!o] Theorems [T} [2] are adapted to a Joo analogue of the situation in Example [2} 
where the x{ti) are specified at various ti and x^^^ is unconstrained. As in Corollary 
[2| a necessary condition for such a feasible curve x to be optimal is that equation 
([2]) should hold along at least one arc, with an additional conditions when the arc is 
terminal. 



3. X'^-Extremals 

A piecewise-C^ field W : [t^^i] — ^ TM along x E F is variational when W(ti) = 
\/i-W\t=t, = for z = 0, 1. For W variational, a variation /i h> x/i G F is given by 
Xh{t) := exp^(^-)(/iM/(t)), and 

{^x,mh=o = w{t). 

For X G F set Sx ■= {s G [to,ti] : || VtX*^^)(t)||s=i = J(x)}. Then Sx is nonempty and 
closed. 

Definition 2. x E F is an Z/°°-extremal when, for any variational field W along x, 
we have {L{W),VtX^^^)t=sw > for some sw G Sx- 

Lemma 1. If x is optimal then x is an -extremal. 



Proof: From the definitions of Xh and L, we find, for any variational field TV, 

^(l|v4'¥)/.=o = 2(LW,V,x(^)). 

Suppose X G -F is not an L°°-extremal. Then, for some variational field W ^ {L{W)^ VtX*^^ 
at all s E Sx- Therefore, and because Sx is compact, for h sufficiently small we 
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have II Vtxjj^^ 11^^^ < J(x)^ for all s G 5*^;. Then || Vta;!^^"* |p < J(a;)^ at all t in an open 
subset U of [to,ti] containing S^. In particular, if U = [to,ti] then x is not optimal. 

If f/ C [to,ti] then, because the nonempty set [to,ti] — U is compact, sup{||Vtx|j^^||^ : 
t ^ U} is \\Vtx[l'^\\j^t^ for some ^ U. Because S.^ C U, \\V tx^^^ < J{xf. So 
||Vjxjj^''|p < J{xY on [to,^i], and x is not optimal. This proves Lemmajl 

Lemma 2. If x is an -extremal then || VtX*-^''||2:(t) is constant. 

Proof: If X is a geodesic the lemma holds trivially. For x not a geodesic we have 
J{x) > 0. Arguing by contradiction, suppose that the open set U := [to,ti] — Sx 
is nonempty. Then U contains an open interval of the form (t^, — 5, t^, + 6) where 
G [to,ti] and 6 > 0. Define fields Wl '■ [toyt* — 5] — )■ TM along x|[to,i* — S] and 
Wr : [t^ + 6,ti] — )■ TM along x\[t^ + 6,ti] by solving the linear initial and terminal 
value problems 

L{Wl) = -VtxW subject to WM = VtWL\t=to = 
HWr) = - Vtx(^) subject to WRiti) = \/tWR\t=t, = 0. 

Then Wl and Wr are the restrictions of a single variational field W : [to;^i] ~^ TM 
along X, for which {L(W),VtX^^'^)t=s = —J{xy for all s G Sx- Since x is an L°°- 
extremal we have a contradiction. This proves Lemma [2] 

Incidentally, there is also a kind of converse. 

Proposition 1. Let x G -F satisfy the necessary condition of Theorem^ Then x is 
an L°° -extremal. 

Proof: Set X(t) := (j){t)\/ tX^^\t) . If x is not an Z/°°-extremal then, for some varia- 
tional field W, 

{VtW + R{W,x^^'^)x^^\X) > on[to,ti). 

Using symmetries of the curvature tensor [13j (2.27), and equation (|2|, the left hand 
side is 

(v?iy,x) + (iy,i?(x,x(^))x«) = (ViW,x) - (iy,v?x) = rW(t) 

where r : [to,^i] -> M is given by r(t) := {VtW,X) - {W,VtX). So r is strictly 
increasing. 

Because W(ti) = \/tW\t=ti = 0, r(tj) = for z = 0, 1. This is a contradiction because 
r is strictly increasing. So x is an L°°-extremal after all, and this proves Proposition 

m 
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Proposition [T] is not used for the proof of Theorem [TJ Instead we need Lemma [3] 
below, which asserts that an optimal curve is also locally optimal. This would be 
trivial if J{x) was the integral of a non-negative function of x and its derivatives, 
as with geodesies, elastic curves and Riemannian cubics. In the present situation a 
proof is needed. 

Lemma 3. Let x G Fxg^xi,vo,vi,to,ti be optimal, and choose sq < si in [to,ti]. For 
i = 0,1 set Ui = x{si) and Wi = x^^\si). Then the restriction y = x\[so,si\ of x to 
[sq,Si] is optimal, considered as a curve in -Fj;o,?/i,'«o,wi,so,^*i ■ 

Proof: By Lemma [2| J{y) = J(x). If y is not optimal let y G -Fyo,i/i,«'o,«'i,so,si be a 
curve with || Vt|/*-^^||s < J{y) for all s E [so,Si]. The track-sum x : [to,ti] — > M of 
x\[to,So], y and x|[si,)f:i] is piecewise-C^ . Indeed x G -Fa;o,xi,Do,t)i,io,ti- 

If So = to and si = ti the Lemma holds trivially. So we can suppose without loss of 
generality that ^2 ^ [sc^i] for some t2 G [to)^i]- So J{x) = \\^tx^^^\t=t2\\x{t2) = J{x). 
So X is optimal. Since || Vfy'-^'^-'lls < J{y) for all s G [so,Si], ||Vtx||;j(f) is nonconstant, 
and this contradicts Lemma [2j So ?/ is optimal after all, and Lemma |3] is proved. 

Given sq = t^ < Si < . . . < Sj < . . . Sn = ti, any curve x G F is a track-sum of its 
restrictions to the subintervals . If a; is optimal then, by Lemma [sj so are 

the with respect to the values and derivatives of x at and Sj. So it 

suffices to prove Theorem [T] when x maps into a coordinate chart of M. In this case, 
M can be replaced by an open subset of M*" equipped with some Riemannian metric 
( , ). From now on suppose this has been done. 



4. Local Geometry 

Mainly to be clear in future on matters of notation, we briefly review some 
coordinate-based differential geometry, where M is taken as an open subset of M™'. 
Readers may prefer to skip to Lemma |4| then ^ and onwards, referring back as 
necessary. 

For any f G M C M'" we have a possibly non- Euclidean inner product ( , on 
j^m _ The associated norm is denoted by || and the dual cD G M™" of a 

linear form u G (M*")* with respect to ( , )„ is given by uj{w) = {Cj,w)v for any 
w G E"". 

Any vector w G M"^ can be written in the form w^Ci, where summation over i is un- 
derstood, where the w'' G M, and where 61,62, ... ,em are the standard basis elements 
of R"^. For any vector field 1^ on M, = W'et where the : M ^ R. The 
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Levi-Civita covariant derivative VyVT of at y G TM^ is the vector in TM„ = 
given by 

dW ■ 

^yW = —Y' + V,iY,W) 

where, for each v G M, : TM^ x TM.^ — )■ TM^ is the symmetric bihnear form 
given by 

V,iY,W) ■= rJ(i;)FW^efc 

and the Christoffel symbols F^^- are functions for all 1 < i,j,k < m. The 
Riemannian metric determines the Christoffel symbols according to formula (2.26) 
of [13], and thereby the Riemannian curvature Rv{X,Y)Z which is trilinear in 
X,Y,Z e TM^ and given by 

= ^ 9^ " 9^ + ^'^^'^^ - Vh)X'Y''Z^e. (3) 
according to formula (2.20) of [13j, where summation over i,j, k,p,q is understood. 

Now, given x E F = Fxf^^xi,vo,vi-i the velocity y := x'^^'> is a vector field along x, and 
the covariant accleration of x is given by the formula 

V,x« := V,a,(,)2/ = y^^\t) + T{x{t)){y{t),y{t)) 

and an optimal curve is one for which the maximum Riemannian norm of the right 
hand side is minimised. By Lemmas [ij [2| it suffices to optimise over curves x E F 
for which z := || ViX*-^-* [[^^{t) is constant. Since Theorem [l] holds trivially when z = 0, 
take z > 0. 

Let S := S"^^^ C E"^ be the unit sphere with respect to the Euclidean norm || \\e 
and define / : M x S M"" by 

II ^ II UI 

Lemma 4. For any u E S and any j = 1, 2, . . . m, 

d 



dwj 



f{w,u) = -{T{ejJ{w,u))J{w,u))^f{w,u). 



Proof of Lemma |4} Because the Levi-Civita covariant derivative is compatible 
with the Riemannian metric, 



did d 
which proves Lemma |4} 
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5. Local Reduction to Optimal Control 



For M open in M*", and any fixed w G M, tlie assignment u i— )• f{w,u) is a 
diffeomorphism from S onto tlie unit spliere witli respect to tlie Riemannian norm 
II So, for any nongeodesic x E F with || Vta;*^"^)(i:)||a;(f) constant, there are unique 
curves y : [to,ti] — )■ M*", u : [to,ii] — ^ S and 2; : [tc^i] — ^ (0, c>o) satisfying 



= y{t) (4) 

= zf{x{t),u{t))-T{xmyit),y{t)) (5) 

= (6) 

with x{to) = xo, x{ti)=xi, and ?/(to) = -^^o, vih) = Vi. (7) 



So an optimal curve x corresponds to a control extremal, namely an extremal {x, y, z) : 
[to, ti] Mx R"^ X (0, 00) with control u : [to, ^i] — ?■ S for the optimal control problem 
with dynamics Q, ([S]), ([6]), boundary conditions ([T]), minimising 

(ti - to)z = [ z{t) dt. 

Jto 

The admissible controls are taken as the piecewise-continuous functions u : [to5^i] ~^ 
5*. This is equivalent to varying x : [to, ti] through piecewise-C^ curves with || ytX*^"*^^ \\x(^t) 
constant, while satisfying x{ti) = Xi and x^^\ti) = Vi for i = 0,1. 

As a control extremal, [x, y, z) and u satisfy the Pontryagin Maximum Principle [25], 
stated in terms of the Hamiltonian H : {0, 1} x (M x x (0, 00)) x (M"' x K"" x 
M)* X 5 -)■ M where 

H{e,x,y,z,X,fi,iy,u) := - ez + X{y) + ij,{zf{x,u) - T^{y,y)). 

The Pontryagin Principle concerns a number e G {0, 1} and a curve of costates 
{X,fi,u) : [to,ti] — )• (M™ X M™ X m)* associated with the control extremal {x,y,z), 
satisfying 

A(^) = fiiz{T{ej,u),u)^u + dT^{ej){y,y))e* (8) 
^(1) = -X + 2fi{T,{e„y))e* (9) 
= e-ij{u) (10) 

for almost all t G [to7^i]- Here we use Lemma |4] to differentiate u := f{x,u), e* is 
the dual of ej, and we sum over j = 1,2, ... ,m. Moreover, if e = then (A, z/) is 
nowhere trivial. 
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Corresponding to the absence of conditions on zito) and z{ti) are transversality 
conditions 

z/(to) = Hti) = 0. (11) 

The Pontryagin Principle asserts that, for almost every t G [to, ti], u{t) is a maximiser 
of H{e, x(t),y(t), z, X(t), fi(t), z/(t), m), namely whenever fi(t) ^ 

/i(t) 

6. Local Proof of Theorem [T] 

We continue to suppose that x is an optimal curve in an open subset M of E™". 
Lemma 5. e = 1 and ||/i(t)|U(t) ^os mean 1. 



Proof of Lemma 5 By (10) and transversality, ^^^^ ||/i(^)IU(t) <ii = e(ti — to). So 
if e = we must have fi identically 0. Then i/ = by (10) and transversality, and 
A = by But for e = 0, (A,/i, i/) is nowhere trivial. So e = 1 after all, and 

is proved. 



l!o IIa^WIUw dt = ti- to. So Lemma 



Lemma 6. L{fi) is identically 0. 

Proof of Lemma [6| For t E [to,ti], rewrite ([S]) and ^ as 

VtA = {fi{z{r{ej,u),u),u + dr,{e,){y,y))-X{Vtej))e* (12) 



-X + fi{Vtej)ej. 



(13) 



Substituting for A from ( 13 ) in ( 12 ) 



V^A = {fx{z{T{e,,u),u),u + dT,{ej){y,y)-el{Vtej)Vtek) + Vtifx)iVte,)e*. (14) 



Differentiating (13) and substituting for VtA from (14), we find 



- V^A + (Vt/i)(Vte,-)e* + /i(V?e,)e* + fi{Vtej)Vte* 
-IJ,{z{T{ej,u),u)^u + dT^{ej){y,y) - V^ejOe*. 



Now dT^{ej){y,y) - V^ej = -^ypy^ek - Vi(rJ.?/pefe) 



p.q ^ pj 

-ypVqek 



dxj 
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OXj OX 



q 

by Because fi{u)u = \\fi\\xU = /i, we have fi{{T{ej,u),u)xU — TpjUpCk) = 

{r{ej,u),u)x\\xp,\\x - fi{r{ej,u)) = {r{ej,u),fi)x\\x - fi{r{ej,u)) = 0. 
So Vffi = -fi{Rpjgypygek)e*, and 

V?/i = - {fi,Rpjqypyqek)xej = -{R{ej,y)y,fi)xej = -{R{fi,y)y,ej)xej = -R{fi,y)y 

using symmetries of Riemannian curvature [I3] (2.27). This proves Lemma [6] 

So yU is a piecewise-C^ field defined along the optimal curve x. Defining $ : [^0,^1] — ^ 
[0,00) by 

m ■■= \\m\i(t), 

we have L{^y^{t)\/tX^^^) = L{zjj) = zL{jl) = 0, since z is constant and by Lemma 

El 

Lemma 7. $^^(0) is a discrete subset of \tQ^ti\. 

Proof of Lemma [7| Let s* be an accumulation point of $^^(0). Then ^(s*) = 
$'•^^(5*) = = (the second derivative is taken to be one-sided if s* is one of 

the finitely many points of discontinuity of <l>*^^^). So 

<P ^ ^ d ^ ^ 

^(/i,/i)x(s,) = 2— (Vt/i,/i)^(,.) = 2(ViV,/^)x(s*) + 2(Vi/i, Vt/i)^(,.) = 0. 

So jl{s^) = (Vt/i)(s=i,) = 0. Since L is a second order linear differential operator, fl is 
identically 0, by Lemma |6j This contradiction of Lemma [5] proves Lemma [7} 

Since fl is piecewise-C^, so is VjX*^^^ on any open interval / C [to)^i] where $ is 
nonzero. So x\I is piecewise-C'*. So fl\I is also piecewise-C^, and so on. In fact x\I 
and are C°°. So ^^'^^ is continuous except where /i is 0. So x and jl are C°° 
except where $ is 0. 

Theorem [1] is now proved for the special case where M is open in M™. As argued at 
the end of ^ this completes the proof of Theorem [T] in general. 
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7. Spheres 

Take M to be the unit m-sphere S"™" in E"^~^^. The Levi-Civita covariant derivative 
VyZ of a vector field Z on S"^ in the direction of F G TS^ = y-^ is given by 

VyZ = dZy{Y) + {Y, Z)y 

where ( , ) is the Euchdean inner product, Z is treated as a map from S"^ to E"^^^, 
and y G S^. Accordingly, the Riemannian curvature R on S"^ is given by 

R{X, Y)Z = (r, Z)X - {X, Z)Y and, for a field X along x : [to, ti] ^ 5'", 

L(X) = X(2) + (x(i),x(i))X + ((x(2),X) + 2(x(i\X«))x. 

If X is optimal then, by Theorem [T| for some z G M and some field X with a discrete 
set of zeroes, we have 

= zX/||X|| - (15) 
X(2) = -(x«,x(i))X-((x(2),X) + 2(a;«,X«))a; (16) 

wherever X(t) 7^ 0. 

Example 4. Ta/ce m = 2, z = 1.2, and 

x(0) = (1,0,0), xW(0) = (0,1,0), X(0) = (0,1,200), X«(0) = (-1, 2, 1). 

The solution x : [0,8] — ?■ S*^ zs shown, viewed from (2,2,-1), in Figure^ beginning 
at the left, and ending at x{8) ^ (-0.433207,0.898726,0.0679917). 

8. Bi-Invariant Lie Groups 

Take M to be a Lie group with a bi-invariant Riemannian metric ( , ) , as in [15] 
§21. The (left) Lie-reduction of a field Y along x : [to5^i] — G is the curve Y^ in the 
Lie algebra ^ of G given by 

y,(t) := dL-l,^Y^^,) 

where L{g) : G ^ G denotes left-multiplication hy g E G. Given xq, xi, fo, fi, and an 
-Fa:o,a;i,t;o,i;i,io,ti"OP^i™^l curvc X : [to; ^i] ~^ C, dcuotc the left Lie reduction of x^^^ by 
V. By Theorem [1) z = \\V'^^\t)\\ is constant, and Xl = (pV^^'^ where X := ^VjX^^). 

Since X/^ is a multiple of comparison of ^ with Lemmas 1, 2 of [IS] gives 

Xf = [Xl'\^] (17) 
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Figure 1: x : [0, 8] -h- 5^ in Example [i] 



namely X^^^* defines a Lax constraint [20] on V. So by Theorem 2.1 of [2T], if G 
is semisimple, x is generically obtainable by quadrature from x|^^ = 0(^^V^(^) + 
01/(2) and V. Also c := ||X[^^||2 = z^(j)W{t)'^ + (f){tf\\V^^\t)f is constant, by 
Integrating (17), 

Xi^) = [Xl,V]+C (18) 
where C G ^ is constant. Taking inner products of both sides of (18) with V^^\ 

z(j)^^) = (C, =^ z(P{t) = {C, V{t)) + a 

where a is constant. Taking inner products of both sides of ( [l8| with V, 

{X^\V) = {C,V) = z<j>-a. 



By analogy with Riemannian cubics the optimal curve x is called null when 
C = 0. Then, since Xl defines another Lax constraint on V, HXj^p = z'^(j){t)'^ is 
constant. So (p is constant unless x is geodesic. Consequently, by equation ^ in 
Theorem [T| we have the 

Corollary 1. For M a bi-invariant Lie group, x G FxQ^xi,vo,vi,to,ti is a null optimal 
curve if and only if x is a null Riemannian cubic. 

Null Riemannian cubics are studied in [TH], [20], |22] for G = S0{3), but most optimal 
curves are non-null and have different geometry, as seen from their reductions. 

Example 5. Take G = 50(3), with the bi-invariant Riemannian metric for which 

ad: ^ so(3) 
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is an isometry, where Euclidean 3-space is viewed as a Lie algebra with respect to 
the cross-product x . Then the Lie reduction of an optimal curve x^^^ can be identified 
with a curve V in E^, satisfying 

l/(i)(t) = (19) 
= WxV + C (20) 

where z ER, C E E^, \\ \\ is the Euclidean norm \\ \\e, and W{t) := ^^(t) E E^. As 
noted previously, there are conserved quantities 

c = llW^Wf, a = z\\W\\-{C,V) = z\\W\\ - {W^^\V). 

Taking z = 1.2, C = -(2, 1, 0), 1/(0) = (1, 2, 3), W{0) = (-1, -4, 6), a numerical so- 
lution of ( 19), (20) for V : [0, 700] — ?■ E^ is shown in Figure^ below. The inner dot is 



at V{0) and the outer is at the other endpomt 1/(700) ^ (2.36765,4.69752,8.40276), 
and the image ofV is viewed from (10,8,0). The curve is nonplanar and seems to be 
approximately periodic. Similarly, cf) = \\W\\ : [0, 700] — M seems to be approximately 



Figure 2: Lie Reduction of x^^) : [0, 700] SO{S) in Example [H] 

periodic, and is apparently everywhere nonzero. The graph over the interval [0, 55] 
is shown in Figure^ If C is replaced by (2,1,0), the reduction V of x^^^ appears to 
be unbounded. The image ofV\[0,5] shown in Figure^ resembles that of a null Lie 
quadratic, spiralling curve inwards towards an asymptotic line. However x is not a 
null Riemannian cubic, because C 0. The lower dot corresponds to V{0) and the 
upper dot is at V{5) ^ (1.77133,4.50895,7.05963). 
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Figure 3: in Example [5] 




Figure 4: A nonperiodic reduction of a non-null optimal curve in Example [5] 
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9. Two Points and One Velocity 

So far, for feasible curves x : [to^^i] M, x{ti) and x^^\ti) have been prescribed 
for i = 0,1. With no conditions on the x^^\ti), J would be minimised by a geodesic. 
Consider next what happens when the condition on x^^\ti) is lifted. 

Definition 3. Let xo,xi G M and vo G TMxq be given. Let Fxg^xi,vo,to,ti be the set 
of piecewise-C^ curves for which x{ti) = Xi for i = 0,1 and x^^\to) = Vq. Defining 
J '■ Fxo,xi,vo,to,ti ^ in the same way as before, a minimis er x G FxQ^xi,vo,to,ti of J is 
said to be Fxo,xuvo,to,ti-'^pi^'^^^- 

Theorem 2. Let x be FxQ,xi,vo,to,ti-optimal. Then \\Vtx^^^\x(t) is constant. For some 
piecewise-C'^ function $ : [to,ti] — )■ [0, oo) with $^"'^(0) a discrete set containing ti, 
we have 

L{(j)VtX^^^) = (21) 
where (f){t) := A/$(t). Moreover x and are C°° except possibly where = 0. 

The change in the conclusion from Theorem [l] is that = 0. We omit the 

changes of detail needed for the proof of Theorem [2j Considering track-inverses, a 
similar result holds when xo,xi & M and vi G TMx^ are given, and feasibility means 
x{ti) = Xi for i = 0,1 with x^^\ti) = vi. 

10. Many Points and No Velocities 

Another kind of feasibility occurs when, for n > 3, we are given real numbers 
to < ^1 < • • • < < • • • < and points xq, xi, . . . ,Xj, . . . ,Xn G M. Denote by 
Fxo,xi,...,xn,to,ti,...,t„ the set of piecewise-C^ curves x : [to,tn] M satisfying x{ti) = xt 
for i = 0,1, . . . ,n. Define J : Fxo,xi,...,x„,to,ti,...,tr, ^ M by 

J(x) := max{||VtX^^^||^.(i) : to < i < in}- 

Then a minimiser x G Fxo,xu-,xr„to,tu-,t„ of J is called Fxo,xu...,xr,,to,tu...,t„-optimal. On 
considering optimality of track-summands, we find from Theorems [T| |2} 

Corollary 2. Let x be FxQ^xi,...,x„,to,ti,...,tn-optimal. Then, for some j = l,2,...,n 

and all t G [tj_i,tj], \\VtX^^^\\x(t) is constant. Moreover, 

• for some piecewise-C"^ function $ : [tj^i,tj\ — )• [0, oo) with $~^(0) discrete, we 
have L(0Vfx(^)) = for t e [tj-i,tj], where (p{t) := ^/^t). 
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• X and are C°° except possibly where = 0, 

• if j = 1 then $(to) = 0, 

• ifj = n then $(t„) = 0. 

Proof of Corollary: Suppose there is no such j. Then, by Theorems [T| [2| for 
every j = 1,2, ... ,n, is suboptimal with respect to the restriction of Joo 

to Fxit,.i),x(tj),xW(tj-i)M'Htj)- This contradicts optimality of x. 

Example 6. For a complete Riemannian manifold M, choose Xq,Xi,X2 € M not 
lying on a geodesic. Choose real numbers to < ti < t2 and let x be Fxg,xi,x2,to,ti,t2- 
optimal. Set e := Joo{,x). Then £ > 0. 

For some t^ > t2 let y : [^2,1^3] M be a geodesic and set X3 := ^(ts). Let y : 
[^2,i^3] M be be any piecewise-C"^ curve with y^'^^t^) = x^^\t2), y{t^) = X3 and 
\\'^ty''^^\\y{t) < £ for all t G {t2,ts). For instance y could be chosen as a perturbation 
of y, or indeed y itself. 

Then the track-sum x : [to, ^3] M of x and y lies in -Fa:o,xi,a:2,x3,to,ti,t2,t3 ■ Moreover 
Jooix) = Joo{x), and x is easily seen to be Fri;Q^xi,x2,x3,to,ti,t2,t3-optimal. Yet, since 
there are continuously many possibilities for x, the conditions of Corollary [1| do not 
apply for j = 3. 

11. Conclusion 

This paper proves necessary conditions for minimum acceleration curves in 
Riemannian manifolds M, analogous to the conditions in [17] for Riemannian cubics. 
Results of this kind were previously known only when the manifold M was fiat. In 
the present paper the necessary conditions are examined in detail when M is a sphere 
or a bi-invariant Lie group. Examples are given in the case where M is S"^ or 5*0(3), 
raising questions about the asymptotics and symmetry of optimal curves. 
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